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ABSTRACT 

An analytical approximation Is generated which follows the 
development of an adiabatic spherical blast wave in a homogeneous ambient 
medium of finite pressure. At early times, when the external pressure Is 
negligible, the structure is that of the usual self-similar solution. At 
later times, the structure evolves smoothly as the shock weakens, the 
post-shock compression declines, and the gradients in pressure and 
density become less severe within the shocked region. The complete 
structure should be reliable down to a post-shock compression of about 2, 
with conditions close inside the shock remaining well-described somewhat 
longer. 

An analytical approximation is also presented for the 
electron-temperature distribution resulting from coulomb collisional 
heating. It is shown that thermal conduction, limited by saturation at 
early times, fades in importance just as coulomb-collisional heating 
becomes significant. An estimate is made of the nonequilibrium cooling 
coefficient and the degree of ionization equilibration expected by the 
time significant cooling sets in. From this, estimates of the endpoint 
of the adiabatic era, based on the collisional equilibrium emissivity, 
are shown to be reasonably accurate. 

The dynamical, thermal, ionization, and spectral structures are 
calculated for blast waves of energy E 0 = 5 x 10^ ergs in a hot 
low-density interstellar environment. A formula is presented for 
estimating the luminosity evolution of such explosions, including the 
effects of nonequilibrium ionization. It is then shown that the B and C 
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bands of the soft X-ray background are reproduced by such a model 
explosion If the ambient density is about 0.004 cm"3, the blast radius is 
roughly 100 pc, and the solar system is located inside the shocked 
region. The age of such an explosion is roughly 10^ years. This result 
is almost independent of whether there is appreciable noncoulomb heating 
of the electrons. The M band count rate for the model is considerably 
less than what is observed, consistent with the suggestion that this 
radiation is Independently produced, possibly in the galactic corona or 
beyond. 

Two regions of the sky, however, which have enhanced X-ray emission, 
and which may be such explosion-heated bubbles viewed from outside, are 
both, like the background, brighter than the models In the M band. The 
0V1 column density measurements also provide difficulties: models with a 
homogeneous ambient medium seem to require an ambient (preshock) 
temperature greater than 3 x 10^ K if the medium extends only slightly 
beyond the shocked region, or greater than 8 x 10 5 K if it fills a 
substantial portion of the Interstellar volume. It is suggested that 
evolution in a pre-existing cavity with a strong density gradient may 
remove both the M band and OVI discrepancies. 

Subject headings: interstellar: matter - shock waves 

X-rays: general 


¥ 


4 


Introduction 

This paper consists of four sections merging material from two 
separable but related investigations. In the first section, an analytic 
approximation is presented for the nonsimilar time evolution of the 
dynamical structure of an adiabatic blast wave generated by a point 
explosion of energy F 0 in a homogeneous ambient medium of density p Q and 
pressure p Q . At early times, this solution resembles the zero pressure 
similarity solution of Taylor (1950) and Sedov (1959), but extends the 
range of investigation well Into the regime in which the external 
pressure is significant. 

In the second section, the thermal history of tiie material is 
discussed. In particular, a scheme is presented for evaluating the 
electron temperature distribution for the evolving structure of §1 If 
only coulomb collisions between electrons and ions serve to heat the 
electrons. Approximations to the distribution at early times are also 
given in order to provide a good description of the rate at which the 
equilibration between electron and ion temperatures takes place. 

The estimate of the cooling time of a blast wave Is reviewed. Such 
estimates, however, are based on the assumptions that the electron and 
ion temperatures have equilibrated and that the Ionization states are in 
co 1 1 1 s 1 ona 1 equilibrium with that temperature. The first of these Is 
shown to be valid before significant cooling occurs, but the second is 
less easily tested. Hence an estimate of the cooling coefficient is made 
for the assumption that the ionization level is lower than in collislonal 
equilibrium, and It Is shown that the cooling rate, although generally 
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higher than In equilibrium, approaches the equilibrium rate before 
significant cooling takes place, verifying the approximate validity of 
cooling time estimates based on the above assumptions. This Is then 
rechecked by estimating the degree of ionization equilibration versus 
time, evaluating the results at the cooling epoch. 

In the third section, a procedure Is presented for following the 
state of a given fluid element trough the evolving dynamical and thermal 
structures of §1 and §11, and a computer program Is described which then 
calculates the Ionization state, emisilvlty, spectrum, luminosity, etc. 
for such a blast wave. The luminosity evakitlon for a set of diagnostic 
runs Is presented and compared with estimates based on §11. In addition, 
for two runs, the dynamical, thermal, and oxygen ionization structure are 
presented, as are the surface brightness spectra. 

In the fourth section, the models generated are applied to a 
discussion of the soft X-ray background. It Is hypothesized that this 
background radiation derives in part from the thermal emission of a very 
large blast wave within which the solar system Is immersed. Using 
equilibrium models and clues from the properties of the background as a 
guide, the required blast wave is estimated to have a radius of about 100 
pc, a characteristic temperature of about 10® K, and to be propagating 
Into an ambient density of about 0.004 cm“3. the required energy is 
about 3 x 10®® ergs, coincidentally similar to the energy released by a 
supernova explosion. 

The diagnostic runs of section § 1 1 1 had been chosen to test this 
hypothesis. Thus the resulting surface brightness spectra are translated 
Into the band count rates which would be observed by the Wisconsin soft 
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X-ray astronomy group's rocket payload. These are then compared with the 

observed count rates, with the result that the two lowest band rates can 

easily be produced, but the medium-energy band apparently cannot. Two 

possible remedies for this deficiency are discussed. 

The OVI column densities calculated for these runs are then 

interpreted and compared with an upper limit imposed by observations. It 
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is found that, in general, there is too much 0 and several ways to 
avoid this are discussed. 

Finally, in §V, the conclusions pertinent to the application in § IV 
are presented. 

Readers interested primarily in fluid mechanics should concentrate 
or, §1 where the extrapolation of an approximation scheme due to Kahn 
(1975) has provided considerable insight into a problem without exact 
analytical solution. Those whose interest focuses on the evolution of 
normal supernova remmants should skim § 1 1 with particular attention to 
I ID . Those who care most about the soft X-ray background, the local 
interstellar medium, or the possibility of a galactic corona could begin 
with §IV. 

I. The Dynamical Evolution of An Adiabatic Blast Wave Into a Region of 
Finite Pressure 

A. Equations of Motion and Normalized Variables 

The evolution of an explosion of energy E Q Into a cold (zero 
pressure) homogeneous medium of density p Q is described by the similarity 

solution of Sedov (1959). The evolution of a similar explosion into a 
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region of non-negllglble ambient pressure, p Q , has been approximated by 
Gaffet (1978), Kahn (1975) has given a method for generating simple 
high quality analytic approximations to the Sedov structure; we here 
generalize this method to give comparable results for the p Q f o case. 

The equations of mass, energy, and momentum conservation In a 
spherically symmetric adiabatic gas with y • 5/3 are, respectively. 
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where R Is the distance from the explosion site, u ■ DR/Dt Is the 
velocity of a parcel of gas at R, and p, p are density and pressure. 

The region of the shock disturbed material is bounded by the shock 
radius, R s , expanding outwardly with velocity v $ . The jump conditions 
relating conditions just Inside the shock to those outside are 



5 ^ 
3p o 
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( 2 ) 


where the normalized variables x ■ p/p q , and y ■ p/p Q have been 
Introduced and the subscript s refers to post-shock conditions. Three 
further normalized variables will also be of use: 


r-R/R,, < • y/(x) 5//3 , U - M(R)/M(R $ ) - / 3r 2 x dr (3) 

where M(R S ) ■ 4 7T R*P 0 /3 and < is the adiabatic invariant for a parcel of 
gas. 
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B. Boundary Conditions at r ^ 0 and r ■v 1 

Early in the evolution of the explosion, the Sedov solution will be 
appropriate, in which case 
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As the evolution proceeds, material shocked early enough to be 
governed by the Sedov equations will be left behind by the shock front, 
eventually becoming part of the central region. In this regie the 
sound speed remains quite high (if the adlabatlclty assumption remains 
valid) and the pressure is nearly uniform. During the flow, any given 
shell of this Interior material retains two constants of the motion, 
the mass enclosed by the shell and the adiabatic invariant, k , the 
latter being determined by equation (4). From the dependence of these 
two constants on one another, the central density structure can be 
Inferred as a function of position and the value, y^ , of the normalized 
pressure In the central plateau. The results can best be expressed in 
terms of a characteristic radius defined by 
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at which the effects of the external pressure begin to be felt. (This 
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particular choice simplifies equat1on(ll ) . ) 

The resulting central density distribution Is 


x - 4(f ? ) 3 ' 2 ,, 3 ' 2 <£ ) S ' 2 

c 


0.30 y, 3 ' 2 (f c )? /2 


(6) 


The particular result from this discussion which Is useful In 
applying Kahn's method to the structure is that, so long as there is a 
central pressure plateau, y is proportional to r^72 making 


J 3y_ , 15 

y 3r 2r 


{central regions). 


(7) 


Oust behind the shock front, the derivatives of pressure, density, 
gas velocity, and y can he found by combining the equations of lotion (1) 
and the jump conditions (2). However, they also depend on the 
deceleration rate of the shock front itself. This deceleration rate can 
be expressed in terms of the parameter 

a = R $ dy s /dR $ (8) 

since knowing a(R $ ) is equivalent, through the jump conditions, to 
knowing v (R ), and therefore R ( t) by normalizing to the Sedov solution 

S S 5 

3 3 

which is appropriate for « R c< In this limit, R & y $ 9R C and a ■* 
-3y s . 

The resulting derivatives evaluated at r ■ 1 are: 
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Abbreviating the partial derivatives with primes, the derivatives at 
early times (when y s » 1, x $ ■ 4, a R -3y $ ) are y' * 8y $ , x' ■ 48, u' « 

2u s , y' = 12, y" » 163. 

C. The Evolving Shock Strength 

In order to proceed further, a reliable approximation Is needed for 
a(R ), or a (y ). We suggest that 

d j 

a * -(y s -l) (3-2/y s ) (10) 

has several desirable properties. First, it approaches -3y s at large y $ , as 

required by the Sedov solution. Secondly, it degenerates to -(y s »l) as y s -*l, 

appropriate to a spherically symmetric sound wave with amplitude Inversely 

proportional to R $ (l.e. R s (y s *l) approaches a constant). Thirdly, aty s »6, 

It has the value -40/3. At y »6, x =2.5, and with this value of a, 

d(lnx)/d(lnr) « 9/2 at x=l as well as at x=0. The density distribution then 

9/2 

has the simple form x=2.5 r at this epoch. Thus the density distribution in 

a log-log plot passes smoothly through a single power law as it changes from 
concave upward at early times to concave downward at late times. 

Finally, since this approximation to a determines the time evolution of R s , 
and underlies much of the remaining discussion of this paper, Eric Jones of 
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Los Alamos was asked to run a purely adiabatic, numerical -hydrodynaml cal 
calculation of an explosion Into a medium of finite pressure.. He did so 
(Jones, 1978), following the explosion to y ■ 2. Because of the finite 
cell size In the numerical calculation, neither the shock position nor post 
shock pressure Is known exactly. Within the uncertainties, however, there Is 
essentially perfect agreement between his results and the integrated form of 
equation (10): 


« s 3 (y s -i) 3 

(3y s -2) Z 


» constant 


2.02 l 

o 



( 11 ) 


giving the evolution of the post-shock pressure with R $ , normalized to the 
Sedov solution at early times. 


D. The Evolving Dynamical Structure 


1. The Kahn Method and its Extension 

Kahn (1975) suggested that since ( 1/u )(9y/3r) 15/2r as r ♦ 0 and 

du/dr ■ 3x at r ■ 1 where \i ■ 1, the general mass distribution might be 
fairly well described by the form 

(l/u)Ciji/3r) - 15/2r ♦ (3x $ - 15/2) r 6 (12) 

where the unknown 8 is found by forcing y" * 3x s [2 + x'/x s J at f * 1 as 
required by equation (9). 

This procedure automatically gives the correct density distribution 
near the edge. Use of the adiabaticity then assures the correct pressure 
distribution there. So long as 8 > -1, it also gives the correct slope 
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for In p near the center, although the density normalization Is not 
guaranteed. In fact, however, the procedure seems to give a quite 
accurate description for a blast wave Into negligible pressure. (See Cox 
and Franco, 1981, for a comparison.) 

We now attempt this procedure for the non-negllglble p Q case. 
Differentiation of equation (12) yields (at r ■ 1) 

y" * - 15/2 ♦ ((3x $ - 15/2)6 ♦ (3x s ) Z (13) 

Setting this equal to 3x s [2 + x ' /x g 3 at r ■ 1 with x'/x $ from equation (9), 
a in that equation from equation (10), and eliminating v„ In favor of x 

using equation (2), 6 can be evaluated. A more convenient parameter for 
later work Is q ■ 8 + 1 and the result Is 


? x (x -l)(7x -13) 
q ' 6 + ’ “7 - -T 4x - S -1) ? 


(14) 


Equation (12) can be Integrated to find the mass distribution 


u(r) * r 15 ^ 2 exp 


[( 


3x s - 


15/2 


) (r q - 1)] 


(15) 


from which the density distribution can be found, using y' « 3r 2 x from 
equation (3). The result is: 

x(r) ■= [-f- + (x s - -|-) r q ] r 9/2 exp[(^-^)(r q -l)] (16) 

(where q ■ q(x s ) from equation (14)’, Similarity, the original location, 

Rj, of a mass element now found at R * R $ r is 

R, - R s [y(r)] V3 5 R s r, 

3 

since y * (R^/R s ) from the initially uniform density distribution. 


( 17 ) 
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The pressure In this parcel just after It was Initially shocked follows 
from combining equation (11) and (17) and solving the resulting cubic 
equation for y s 1 (u* y $ ) to give 



where Q(v»,y s ) ■ i.(r) (3y s -2) 2 /(y s -l ) 3 and y $ characteHtel the present epoch. 
The density of the parcel just after being shocked was 

. 4 *s1 + 1 09) 

from equation (2). The adlabatlclty of the flow can now he used to find the 
present pressure and temperature In the parcel: 



T - T Q [y(r)/x(r)3 (20) 

where y^ and x^ are both functions of r, and T Q Is the ambient external 
temperature. The entire structure is thus obtained analytically. An example 
Is presented In Sill. 

A difficulty arises, however, when one wishes to follow a collection of 
gas parcels In time to obtain their Ionization structure. This Involves 
Inverting equation (15) to find the present location cf a chosen mass point. 
Our procedure for doing this Is described in f 1 1 1 A . 
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The time, t(R s ), is found by integrating the differential equation 
between v and R . This equation results from eliminating y from 

S S 5 

equation (11) in favor of v , using equation (2), to yield 


z = 


4 l/3 (15H s 2 ■ 11) 2/3 
5 (M 2 - 1) 


( 21 ) 


where M - v /c and c is the ambient adiabatic sound speed. This 
s s o o 

equation was approximately inverted by curve fitting (to an accuracy of 
better than Q.1% for M ^2) whereupon 



z 3/2 


- 0.0933 z 5 e~ z 


( 22 ) 


was the result. At = 1.5, z = 2.039, the fit gives = 1.519. At M s 
= 1.25, z = 3.030, it gives = 1.224 so that reasonable accuracy is 

maintained even to R s ^ 3R C . The time is therefore 


t(R s> * £ igzT ' t. - °- 0933 25 dz 


z .3/2 


5 -z- 


(23) 


R r o 7 5/2 

- C r C Z 

~ 7" L 5 2.68 
'o 


-0.0933 J(z)] 


where 


J(z) = 120 - (z 5 + 5z 4 + 20z 3 + 60z 2 + 120z + 120) e’ Z , (24) 

completing our approximate but analytical description of the dynamical 
variables. 

The general character and quality of the results are summarized by: 

1) y s » 1. The solutions approach Kahn's very accurate approximation 
to the Sedov solution. 

2) y s decreasing toward 6. The structure Is well behaved and becoming 
less severe. The central pressure plateau is becoming flatter. 
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3) y, - 6. x » 2.5, R /R - 1.27, T/T « 2.4, q • 1.25. At this 

5 5 5 C 5 v 

epoch, p ■ r^ 2 , x ■ 2.5 r 972 , and the central pressure distribution 
Is extremely flat with yj <* 2.01. 

4) 6 > y s > 3. The approximation begins to have difficulties near 

r » o, where a presumably non-physical rise appears in the pressure. 
The description should remain reliable, however, for the vast 
majority of the material. 

5) y. ■ 3, x ■ 1.857, R /R B 1.83, T/T c 1.615, q ■ 0, At this epoch, 

b •> 5 V J W 

U » r 39/7 , X - (13/7)r 18/7 , y + 0.41/r 9/7 as r ■+ o. At this point the 
central pressure rise has diverged to the forlh glfen. This behavior 
arises from the inability of the Kahn approximation schema to follow 
c density distribution which is intrinsically concave downward in a 
log-log plot. Even at this epoch, however, 94 per cent of the mass 
is exterior to the pressure minimum at r = 0.6. In addition, that 
minimum has y * 1.0, making this the earliest epoch at which the 
internal pressure is anywhere below the external pressure. This Is 
confirmed in Jones' (1978) numerical results. Furthermore, detailed 
comparison shows that the structure remains well -described outward 
from about r ■ 0.75. 

6) y «2, x ■ 1.5, R /R » 2.52. This Is the latest epoch at which the 

' 's s s c 

post shock pressure evolution of equation (11) was tested and there- 
fore may be the latest epoch at which the post shock derivatives are 
reliable. At least until this point, the outermost structure should 
be reasonably reliable. The time evolution of equation (23) remains 
satisfactory to this epoch as well. 

In his approach to this problem Gaffet (1978) also used the mass coordin- 
ates to follow the adiabatic material, and derived results equivalent to 
equation (9) relating shock strength, deceleration, and post shock derivatives. 
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Similarly he Incorporated a deceleration rate equivalent to a —40/3 at 
y s ■ 6. The greatest Intrinsic difference In his method Is that he assumed 
that dy/dv* was approximately constant rather than o£n w / d£n r being given 
by a sum of two power laws as In the Kahn method. His approach has the 
inherent disadvantage that the results could not be expressed analytically. 

By graphing y versus u for a broad range of epochs from the results of 
this paper, we find that dy/dv Is In fact approximately constant for the 
outer 95% of the mass, but that the average slope can differ by several 
percent from the Immediate post shock value determined by the acceleration 
parameter. This should lead to modest discrepancies between the two approxi- 
mations, but these could be made to vanish in Gaffet's higher order approxima- 
tion which allowed quadratic terms in dy/dv- 

A second intrinsic difference, however, is the method by which the shock 
deceleration parameter was deduced. In our work, a(y $ ) was deduced ab initio 
from plausibility arguments and checked by a numerical hydrodynamic calculation. 
In Gaffet's simplest approximation, he seems to have been led by energy argu- 

3 

ments to a deceleration equivalent to a »-3y + 14/3, making (y r - 14/9) R 
a constant. This Is very similar to our result for y $ > 6, but clearly 
incorrect at large radii. Since Gaffet presents the post shock temperature 
and density structure evolution graphically for epochs as late as y $ * 14/9, 
an Improved deceleration must have been incorporated. The way in which this was 
done, however, is unclear to us. 

When our results are graphed In forms comparable to Gaffet's Figs. 2 and 3, 

we find that our density structures are essentially Identical as expected for 

those epochs with y > 6. For lower values of y there is total disagreement. 

s s 

At small radii, our results are becoming inappropriate as described above. 
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whereas Gaffet's approach does not suffer quite the same fault. At large 
radii, we would expect continued agreement If the same deceleration rate were 
used. The most extreme example presented has T s /T q ■ 1.2, y s ■ 1.56. Our 
post shock density slope Is approximately twice that of Gaffet. Ours corresponds 
to a ■ -0.96, his to a <* -0.76. These two do not correspond to the same shock 
radius since ours exhibits greater deceleration, so that the two density 
distributions at the same shock radius would be less different, but would have 
different shock strengths. In short, whatever method Gaffet used to find the 
deceleration rate did not quite agree with ours at late times. 

Our results for v $ versus R $ (Fig. t In fcaffet) ard alsl In disagreement. 

The biggest part of this discrepancy derives from the fact that Gaffet 
normalized his results to only a rough approximation to the Sedov solution at 
early times, nrklng the slopes of both his lines In his Figure 2 about 9% low. 

Finally we would like to point out that neither method should be considered 
seriously for a description at late times. As early as y s B 6, a ripple can 
be seen developing In the temperature structure. At very late times, this 
ripple near the edge should be found In temperature, pressure, and density. As 
It propagates outward. It leaves behind a hot cavity (described below) within 
which the pressure Is nearly uniform. Neither method Is capable of adopting 
this form. 

2. The Theoretical End Result 

If Ihe blast wave remained adiabatic indefinitely, the conditions In the 
vicinity of the explosion site long after the occurrence can be found. The 
procedure is similar to that used to find equation (6). One need only assume 
that the distribution of adiabatic constant, <, versus enclosed mass, as 
Implied by equation (11), Is preserved as the system returns to the original 



18 


ambient pressure. The result of this calculation Is shown as x_ ( R/R_ ) In 

f C 

Figure 1. The figure also shows the amount of mass which the blast wave 
removed from the volume bounded by R/R £ . 

From the figure, the net result of the explosion Is seen to be the 
generation of a hot (T^ ■ T Q /x f ) cavity with a steep density profile, leveling 
off around R * 2R £ . About 90 percent of the mass Initially with R < R c has 
been expelled. The amount of expelled mass continues to rise, however, 

3 

reaching about 2.5 M (where M * 4irp R /3) by R * 2.5 R , the latest that 
equation (11) has been tested. This continued rise of thg expelled mass would 
be incorrect If the system did not return identically to the original pressure 
In a finite time. 
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II. The Thermal Evolution 

A. Electron Heating and Thermal Conduction 

A long-standing uncertainty In the application of blast wave theory 
to supernova remnants Is whether the electrons can be expected to have 
the same temperature as the Ions (see McKee and Hollenbach, 1980). If 
only coulomb collisions act to equalize these temperatures, then the 
electrons are much colder than the Ions just inside the shock front, and 
require a considerable length of tim* to heat up. Several authors (e.g., 
Itoh, 1978) have shown that considering the remnant as a whole, the 
electrons would reach equilibirum with the ion temperature over an 
appreciable fraction of the interior mass only after a time 

t eq ■ 50°0 yrE 5 , 3/14 /n 0 4/7 (25) 

where Ejj is the blast energy in units of 10^ ergs. Alternatively 
various plasma Instabilities may cause the electrons to equilibrate, at 
least partially, with the ions in the shock front. We have no Immediate 
preference for one of these two ideas, hoping that the answer will come 
eventually from careful interpretation of x-ray observations of supernova 
remnants. (Again see McKee and Hollenbach, 1980, for a discussion of 
progress to date.) 

A second complication which has drawn considerable attention Is the 
role of thermal conduction in supernova remnants (e.g., Solinger, 
Rappaport, and Buff, 1975; Chevalier 1975, 1977; Cowie, 1977). This is 
particularly important in that thermal conduction, if it were as 
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Important as straightforward estimates Indicate that It could be, would 
ruin the assumption of adlabatlclty which underlies the entire treatment 
of II. 

Such an estimate of the Importance of thermal conduction can be made 

by assuming the adiabatic structure, calculating the thermal conduction flux 

expected, and assessing Its effect on the jump conditions of the shock. The 

modified jump conditions are given by x $ [1 + (2F/p q v s )] * 5 x s -4 where F 

Is the arriving thermal conduction flux per unit area. Thus x ■ 2 follows 
3 3 

from F ■ p v /4, and x * 3 from F * p v / 9 . Clearly a very modest conduction 
os s OS' 

flux is sufficient to alter details of the structure considerably. For the 

8 12 -4 

early adiabatic blast wave, p « R , n « R near the edge, so T « R and 
-3T/3R - 4 T s /R $ . 

Taking k s 6 x 10"^ T^cgs, and defining 6 * 9F/(p o v $ ^) as a measure of 
the significance of thermal conduction, assuming the Sedov evolution, we have 


6 = 


36 K T s 

V,\ 


3.2 


nkv R 
o s s 


0.16 


*51 


3/5 


n 6/5 1 1475“ s 
n o % 


-7 


(26) 


4 

where t^ * t/10 yr. Three conclusions are evident: (1) If 6>>1 In this 

formula, the thermal conduction Is saturated so that the true ratio Is of 
order unity; (2) If 5«1, thermal conduction has little effect on the 
expansion, although Its effects In the remnant interior are somewhat more 
persistent; (3) the changeover from large to small 6, marking the fading of 
thermal conduction significance, occurs precl' 5 tously (measured particularly 
by R s dependency) at a time 


Snech " 5200 * r E 51 3/1 V /? 



when 6 ■ 1 and x $ * 3. 


(27) 
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There are several distinct situations which warrant discussion: 

1. The magnetic field within a remnant may Interfere with thermal 
conduction to a degree that the evolution Is essentially adiabatic. 

2. Electrons may be heated only by coulomb collisions. In which case 
they gain access to the energies of the ions (at t ) only as their 
ability to redistribute that energy Is waning (at t ^). Their 
Influence on the structural evolution Is then never very great, but 
prior to t eC j the ion gas Is nearly adiabatic, the electron gas 
nearly Isothermal, and T fi < T^. After t , the temperature structures 
merge just as adlabatlclty becomes a reasonable approximation for 

the single fluid. 

3. If electrons were heated to the ion temperature at the shock front, 
even for t < t , saturated conduction would lead to a modest 
enhancement of the flux arriving at the front, decreasing x 

and increasing the shock speed somewhat. The greater significance of 
the conductivity In the center would lead to a significant flattening 
of the electron temperature distribution. Once again, however, the 
electrons In the interior would not have access to the ion energies 
for t < t . Hence the flattening of T fi would not lead to a similar 
flattening of either the ion temperature or the density distribution 
(Cowle, 1977). Once again for t < t the ion fluid would be nearly 
adiabatic, the electron fluid nearly isothermal, but the two would 
have the same temperature at the shock. As t t , the electrons 
would be heated by the ions, but the importance of thermal conduction 
thereafter would be small. 
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4. If electrons were not only heated at the shock but some process 
existed to keep the electron and ion temperatures equal thereafter, then 
for t < t both electron and Ion temperatures would be flattened by 
conductivity and there would be a density plateau in the Interior (c.f. 
Chevalier, 1975). Thus at early times, a description based on the 
assumption of adiabatlcity would be inappropriate, although the radiative 
properties would not be drastically different (e.g., Solinger, Rappaport, 
and Buff, 1975; Gronenschild and Mewe, 1979). For t > t mech , an 
adiabatic description would gradually become appropriate. 

In the remainder of this paper, we will consider only two models, 

one with only coulomb heating of electrons, and one with T e * * T 

everywhere. Both models are adiabatic. The former is appropriate for 

possibility 2 above, whether or not the magnetic field quenches thermal 

conductivity. It also approximates possibility 3 for t > 0.6 t because, 

eq 

as shown by Itoh (1978) and SIIC below, this case leads to a flat electron 
distribution which rises to the shock temperature at 0.6 t , providing a 
structure similar to that for possibility 3 for this and later times. 


The second model is appropriate for either possibility 3 or 4 if 
thermal conduction is quenched by the magnetic field, but neither model 
describes either 3 for t < 0.6 t or 4 for t < t^. The two models 
converge, as do all the possibilities, for t > t . 
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B. Itoh's Formulation of Electron Heating 

Itoh (1978) has shown that In an adiabatic flow with T /m » T./m. 

e e 11 

and approximately constant value of In A - ln[1.2 x 10 5 (T^ 2 T e /n^ 2 )] t 

the evolution of normalized electron and ion temperatures, using 

Spltzer's (1962) formula for energy exchange by coulomb collisions. Is 

Independent of the history of a parcel of gas, once Its conditions have 

been established at the shock front passage. Defining T as the average 

particle temperature (as found from adiabatic models which assume equal 

electron and Ion temperatures), Itoh's result is that 

g = T q /T is determined by 
6 


[3 ln ( LWs_ 

2 i - sr 


) - ^T(g + 3)L * in_A_ 
g o 81 




(t 




(28) 


where t Q Is the time the particle was shocked, g Q Is g just behind the 
shock, and n/l3/2 t s * constant established at the shock. The numerical 
form of this equation assumes cgs units, n^/r^ * 0.1. 

In order to find g(t), the above equation must be Inverted. We have 
done this, approximately, by curve fitting with the result that the 
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equation 


f - 4- In (L±- ^) - Sq- (g «■ 3) 
2 1 - v§" 


( 29 ) 


Inverts to 


9 ; 1 - exp [ - (r-f- f) 0 ' 4 (1+0.3 (-§- f) 0 ' 6 }] (30) 

with an accuracy of bitter than 2% f or 0 < g < 1 (although 1-g Is only 
this accurate for g < 0.8). The electron temperature Is thus followed In 
time by evaluating f Q from equation (29), with g ■ g Q , adding the right 
hand side of (28) to f 0 and suostltutlng the result In (30). The 
validity of the procedure depends almost exclusively on adlabatlclty. 

C. Application of Itoh's Formulation to Kahn's Results 

The previous results can be applied directly to obtain the electron 
temperature distribution In the central regions of a strong adiabatic 
blast wave. For this material, the initial Ion temperature was extremely 
high so we take g Q ■ 0. In addition, this material was shocked at a time 
t 0 (R) very much less than the present age of the explosion (In order now 
to be found near the center). 

Hence 


f In A 

center * "51 


T 


TV * 1 


t. 


( 31 ) 
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Also near the center, n/T^/2 is very small so equation { 30 ) can be 
approximated (for f < 0.1) by 

Center = (“T ^ <32) 

whereupon 

~f 'center = I^r • t33 > 

From equations (14) and (1C with x s * 4, the density distribution near 
the center is 



center - 


~T 


-9/16 


, R t9/2 

( nr~ ) 


(34) 


from which the initial location, , of a parcel now found at R can be 
evaluated to have been 



r 3/ ’ 6 ( /-) 5/2 . 


(35) 


The initial mean post-shock temperature of this material was 


T 


1 



T, ,9/16 ( _S,lS/2 


(36) 


where T s is the present post-shock mean temperature. This material has, 
however, been adlabatlcal ly decompressed from 4n Q to n so that its 
present temperature is 
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t - Ti (^-) 2/3 * y - f -> 2/3 . 3/,e ( V /z - 


( 3 ?) 


Substituting equations (34) and (37) Into (33), 


In A _ T 5 5/3 -3/8 * n 2/5 


T e-[-4T n o*s 8 


t r 


(38) 


which result Is Independent of Ri The existence of this central plateau 
In T e has in fact been found previously In the (corrected) numerical 
study of Itoh (1978) and in that of Gronenschlld and Mewe (1979). 
Equation (38), however, predicts the value of this plateau temperature 
analytically through a combination of Itoh's and Kahn's approximation 
schemes. 

Instructive manipulations of equation (38) can be made by 
eliminating either T s or t using the blast wave evolution of equations 
(4). For In A = 35, the second of these yields 


e 'plateau 


* 1.4 x 10 


7 


K (T $ /10 


6 k) 1/15 c 2/15„ 4/15 


-51 


(39) 


showing an extremely weak dependence on T . 

T . t -tjn 

e 


-6/5 

Eliminating T s * t ' shows 


Since T Is dropping more steeply in time than T , there comes a 
time when the nearly uniform electron temperature equals 1 Shortly 
thereafter, a significant fraction of the mass obtains T g = so that 
this moment marks the onset of equilibration. With suitable 
substitutions (as above) in equation (38), 
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T . t r \ Ha t 28/25 

'‘plateau = s [?60 ° * r (40) 

(again with In A * 35). As noted previously, the plateau T g reaches T $ 
at t * 0.6 t . Finally, we note from equation (39) that the electron 
temperature for t < t (and therefore the highest expected at any time) 
depends almost exclusively on the density of the ambler' (assumed 
homogeneous) medium. 

D. Cooling of an Adiabatic Blast Wave 

The cooling rate in and cooling time oi ; a blast wave is a somewhat 
strange thing to contemplate estimating when the electrons are out of 
equilibrium with the ion temperature, and the ionization level 
distribution is even further out of equilibrium with T g . As this section 
will show, however, an equilibrium calculation is reasonable when 
estimating the cooling time marking the end of the adiabatic era. Thus 
we begin with such a calculation. 

-3 -1 

The emissivity of a gas is usually defined as L(T g )n erg cm s 

and values of l(T g ) for a gas whose ionization level is in collisional 

equilibrium with T ore fairly well known (e.g. Raymond, Cox, and Smith, 

1976, hereafter RCS). Kahn ( 1 97 5 , 1976) has pointed out for the special 
- 1/2 

case with L » T g and T ? * T, that a freely cooling parcel of gas has a 
cooling time which is independent of the history.. More specifically, 
p/p reaches zero at a time which is independent of the separate 
histories of p and p . For n * 1.1 and Hp = 1.1 n, this cooling time 

iS 
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At , 
cool 


2.1 kT 

RT)n ’ 


(40.5) 


Furthermore, for the temperature range from 4 x 10^ K to 5 x 10^ K, the 
equilibrium cooling function of RCS is fairly well approximated by L(T) ^ 
1.3 x 10~ 22 (T/10 6 K) ^ 2 erg cm 3 s" 1 so that equation (40.5) gives a 
reasonable estimate of the cooling time for this material. 

Assuming the appropriateness of equilibrium cooling, equation (40.5) 
can then be used (Kahn; 1975, 1976) to calculate the time at which each 
parcel of gas swept over by a stronr blast wave in the adiabatic phase 
will eventually cool. The earliest cooling of any parcel then marks the 
effective end to the adiabatic phase. For a parcel encountered at time 
t s when the post shock temperature was T s , the cooling is complete at a 
time 


t - t + 2J kTs 
cool *s + T* L(T^ 


(41) 


2.1 kT sl 


h <-tr> + <-tr> 


s \ -9/ 5 


where t] is a fiducial time which we choose to equal the shock time of 
the first parcel to finish cooling. Setting d(t coo i)/dt = 0 and t $ = t-j 
at this moment, we find the condition 


2.1 


kT 


si 


5 -T- nrTTO 


= i 


( 42 ) 


and t COQ ^ = ( 14/9 ) t-j . Using the approximate RCS cooling function. 
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h * TT t coo1 = 3 x 10 4 yr ^ 

n o 


T s (t 1 ) = 0.96 x 10 6 K (E 51 n o 2 ) 1/7 . 


Equation (44) verifies that the cooling is occurring at temperatures for 

-3 

which the approximation to L(T) is reasonably valid as long as 2 x 10 < 

2 11 

n Q < 8 x 10 . Finally, comparing equations (43) and (25 ) » * T is 

a very acceptable assumption well before cooling becomes appropriate, so 
that our attention can be concentrated on t’.e effects of the 
non-equilibrium ionization structure. 

We now proceed to an estimate of the non-equilibrium cooling 
coefficient behind a shock front. We note first that the radiation from 
ions of a given element will have a softer spectrum than in the 
equilibrium Cc.se because lower stages of ionization are present. 

Secondly, the net emissivity will be higher than in equilibrium, because 
these lower stages are more easily excited by electron collisions. 

As an element enters the shock front, ic is rapidly ionized through 
the lower stages and soon finds itself caught in a stage for which the 
ionization time is very long. This occurs because the ionization 
coefficients decrease very rapidly with increasing ionization stage. 

Thus the most prevalent stage of ionization will be that for which the 
ionization coefficient is of order l/(n Q t). These coefficients can be 
written (e.g. Cox & Tucker, 1969) 


1.3 x 10' 8 T 1/2 ^ F § % cm 3 s" 1 


<W 


(45) 


30 


where l and F are factors of order unity. From this we can solve for the 
exponential factor 


s - ,/kT e 


__ V 

1.3 x 10' 8 T ' /2 CFn Q t 


(«) 


The emlsslvlty due to that element comes primarily from the colllslonal 
excitation of this Ion at a rate (e.g, Cox and Tucker, 1969) 


E a E ; 2.7 x 


10- 15 T. 


• 1/2 


f 9 .' E/kT e 


erg cirl^ s 


(47) 


where again f and g are of order unity. Thus the cooling coefficient for 
this element can be written 


MV 


-15 -l /? ( lev) 

», (!.. . .. r, 

o c 


iE/I 


(48) 


where Is the abundance of the element relative to hydrogen. The ratio 

of strongest excitation line energy, E, to ionization energy, I, shows 

considerable variation among ions but generally lies In the range 0.25 to 

0.75 and we take E/I ^ 0.5 as representative. We similarly take 0. i to 

1/2 

be characteristic of the quantity fg/(£F) , making 


MV 


^ 2.4 x 10" 12 A, 


I 


ev 

3/4 


O 


1/2 


3 -1 

erg cm s 


(49) 


Finally, realizing that the reason the element became stuck in this stage 
of ionization is that I exceeded kT e , we take I ^ 1.5 kT e (which we can 
safely do when it is not in an exponential). Our final estimate is then 
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(summing over all heavy elements and taking EA^ ^ 2 x 10 ’ 3 ) 

L(T e ) v 0.6 x 10~ 13 T fi 1/4 (ppf) 172 erg cm 3 s" 1 . (50) 

o 

We now eliminate t using equation (4) and, assuming T e ■ T $ , obtain the 
approximate non -equi librium cooling coefficient for material just inside 
a strong adiabatic blast wave with hot electrons; 


L(T S ) * 2 


10" 27 T 


2/3 


-j ) 176 erg cm 3 s ’ 1 


■51 


(51) 


2 

We note that the important scaling of collisional processes with E Q n o 
(c.f. Chevalier, 1974) has survived this estimation procedure. This 
scaling also appeared in equation (44). 

The ratio of equation (51) to the approximated equilibrium 
coefficient is proportional to T $ ^and equation ( 5 ]) is projected to 
drop to the equilibrium value when T $ « T^, t * tg with 

T e n, 5 x 10 6 K (E 51 n Q 2 ) 1/7 

E 3/4 (52) 

t E * 7600 yr 


Checking back on the validity of the approximations leading to this 
result, we find that at time , the quantity in square brackets in 
equation (48) is of order 0.1 so that the error introduced by the 
uncertainty in E/I is no greater than the general factor of 3 or so 
uncertainty from the choice of parameters. The combination of parameters 
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introduced above normalizes equation (50) to agree very well with the 
detailed machine results of Sill. 

Since the true non-equilibrium cooling rate Is always higher than 
that for equilibrium, we can Interpret the above results to Imply that 
for t > t^, the equilibrium rate should be reasonably appropriate. And 
since t^ n, 1/4 tj , we would conclude from this argument that equation 
(43) can be used to reliably predict the end of the adiabatic phase of a 
p » 0 blast wave. 

To test this conclusion, we now approach the question of Ionization 
equilibration more directly. If, at temperature T s , a particular element 
In colllsional equilibrium has z as its dominant ionization stage, and 
that stage has ionization coefficient ctj(z), then this element should be 
close to ionization equilibrium only if 4.4 n Q aj(z)t is large compared to 
1. The dominant stage at any temperature is the one with most nearly 
equal ionization and recombination coefficients. If radiative 
recombination dominates, it is straightforward to show that this would 
imply that I(z)/kT would be about 10 (varying slowly from about 13 at z » 
1 , to 9 at z *10), and that both the ionization and recombination 
coefficients would be roughly 2 x 10"^ *2y-l/2 cm 3 s -l t p or m0$t stages 

p 

of ionization, I(z) < (z+1) 2 I H , where 1^ * 13.6 eV, making (z+1) > 

-5 

I ( z)/I^j ^ 10kT/I H 6.3 x 10 T. Ignoring the difference between z and 
z+1, the characteristic ionization and recombination rates would be 
a (T) -v 1.2 x 10 ^^cm\ . This, however, is a lower limit, owing 

vC 

to our neglect of dielectronic recombination which pushes the equilibrium 
ionization level downward at a given temperature. This sometimes 
considerably enhances the characteristic ionization rates, and shortens 
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the time required to establish equilibrium. We thus take the above value 
as the minimum rate at which equilibrium is established. 

The minimum degrees of equilibration at times t^ and t^ are thus 
measured by 4.4 n Q a^ E (T £ )t E ~ 3(^ 5 i n 0 2 ) 3 ^ 14 and 4 * 4 n Q a^ E (T s (t 1 ) ) * 
5(Egi n o^) 3 ^ respectively. The implication is that at least for 

O 

(^51 n o ' ~ co ^i s ^ ona l equilibrium is nearly valid by the time 
significant cooling begins, supporting the collection of cooling time 
estimates which have been made for normal supernova remnants on the basis 
of equations like (43). 

_ o 

For ambient densities substantially lower than 1 cm , we have not 

so far succeeded in assessing the degree of collisional equilibration 

which occurs before time t-j . We note, however, from equation (44) that 

the post shock temperature at time t-j decreases slowly with decreasing 
2 6 

(E 5 i n 0 ), and is less than 10 K at the densities whose equilibration is 

in question. It is further true that the difficulty with equilibration 

is generally most severe when ionization of the K-shell of abundant 

elements is involved, because of the large jump of ionization potential 

at the helium-like stage. The production of the helium-like stage, 

however, occurs readily owing to the ease with which the lithium-like 

stage is stripped of its last L-shell electron. Finally, for 

temperatures below about 1.5 x 10® K, 10® K, and 6 x 10® K, the abundant 

elements oxygen, nitrogen, and carbon, respectively, are no longer 

ionized beyond the helium-like stage in collisional equilibrium. Hence, 

for T s < 10® K there is a rapid increase in the ease with which 

equilibrium is approached. This should remain true down to E51 n Q 2 'v 1.3 
-5 -3 -3 

x 10 , rig % 5 x 10 cm at which point the complete ionization of 
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helium can no longer be assured; but at such low densities, the effects 
of the ambient pressure and pre-ionization level can no longer be 
neglected. Thus for all cases in which w« would be interested In 
applying the cooling time formula, collisional equilibrium should be 
fairly well approached before significant cooling sets in. Of course 
there will always be a region just inside the shock front in which the 
ions are scrambling towards equilibrium and where the emissivlty is high, 
but the added cooling rate from this layer dwindles with time. 

Finally, by comparing equations (25), (27), (43), and (52) we find 
that, for a strong {p Q * 0) blast wave, t^.tj, t £ , t mech and t are 
always found in the ratios 9.3: 6: 1.5: 1.0: 1.0 as an incontrovertible 
sequence and conclude that equation (43) solved for t coo ^ reliably 
denotes the end of the adiabatic evolution and the onset of shell 


formation. 
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III. Ionization Structure and Luminosity 

A. Procedure for Finding the Ionization Structure 

The previous sections provide good analytical approximations for the 
dynamical structure of an evolving adiabatic blast wave for R < 1.7 R . 
In order to find the spectrum of such an explosion, however, the 
ionization structure must be known at each location and this can only be 
found by following the individual gas parcels' histories from the time 
they are shocked until the moment when the spectrum Is to be evaluated. 
For each such history, the time derivatives of each stage of ionization 
of each atomic species must be integrated numerically. Such computations 
have recently been carried out for strong blast waves (negligible 
external pressure) by Itoh (1979) and by Gronenschlld and Mewe (1979). 

The procedure used here for the integrations follows that outlined in Cox 
(1972) and will not be discussed in detail. The ionization and 
recombination rates used are those of Raymond and Smith (1977). 

We found it useful to choose the parcels to follow, not by the 
radius at which they were located, but by the post shock pressure, y s ^ , 
when they were first struck by the shock. The suitable distribution of 
these pressures was chosen to be sharply peaked near the final post shock 
pressure, y $ ^, at the moment the spectrum is to be evaluated. This was 
necessary in order to resolve the rapidly changing conditions near the 
edge. Thirty parcels in all were followed. 

Once y ^ is chosen, the initial compression, follows from 
equation (19) and the initial location of the mass, (the shock 
position when , » y ^), given by equation (11). The time, t ,, is 

5 51 ^ i 
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found from equation (23). The evolution of each parcel In time was 

followed by subdividing the pressure drop y s ^ - y s ^ into 24 Intervals 

(with smaller Intervals at first and larger ones as y s -► y s f) and 

following conditions In the parcel as the shock moved outward and the 

post shock pressure dropped through the sequence of Intervals. 

At an epoch when the post shock pressure is given by y s , conditions 

in the parcel of interest can be found only after the parcel is located. 

This location follows from knowing R s from equation (11) and since p ■ 

3 3 

(Rj/R $ ) » r, for the parcel, equation (15) becomes 

5/ , (x, - 5/2) 

r 1 « r 5/2 exp [ (r q - 1)]. (53) 


In this equation, and q follow from equations (2) and (14) with the 
current value of y s . The inversion of equation (53) to find 
r(r i ) is facilitated by rewriting it as 


(-f“ y s ' exp ( T” x s ' ^ r l 2q/5 

* C(-f-) x s - exp «-f x s ' l > ^ 


(54) 


a 

which is of the form B a Ae where B is known and A desired. For 4 >_ 
x $ > 2 and 1 > r^ > 0, the range of B is 1.1 > B > -0.2, The inversion 
proceeds easily by iteration with A Q * B/(l + 0.763 B), A 1+l ' 

B exp [-A^ ] . 

In calculating the spectrum, this procedure is also used to locate 
the boundaries of the gas parcels and differences taken to find their 
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volumes. For parcels very close to the shock front, these volumes are 
very small (by choice) so that the number of iterations required to find 
r to sufficient accuracy can be fairly large. 

Once A, and therefore r are found for the parcel, the density Is 
found from equation (16), the pressure and mean temperature from equation 
(20), and the time from equation (23), In models for which the electron 
and Ion temperatures were assumed equal, nothing more was required to 
follow R(t), n(t), and T e (t) for the parcel, and the ionization equations 
can be integrated. We assumed tnat the ions outside the shock were in 
collisional equilibrium with To, the external ambient temperature, as the 
initial conditions for the integration. 

For those models in which electrons were assumed to be heated only 
by coulomb collisions, the post-shock electron temperature was taken to 
be T 0 , so that g 0 » T 0 /T $ .j. On the ^ uture » m P 1an to use a P ost shock 
electron temperature equal to to allow for the adiabatic 

compression.) Then f Q is evaluated from equation (29). At all 
subsequent times. 


f o* 


Ip a 

81 


( t ? 7T ) {t " t si ) 


f + 0.43 
o 


0 / si i V2 

7172 

'0 y si 


t ;i^ 


(55) 


with 1 , A 3 35 appropriate for the casts we explore here. Thus g * T g /T 
is found from equation (30); T is found from equation (20) as before, and 
the determination of T g (t) is completed. The integration of the 
ionization equations can then be carried out. 
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B. Program Capabilities 

The procedures outlined above were written into a computer program 
designed to start with the ambient medium In colllslonal equilibrium with 
temperature T , density n Q , and given elemental abundances, and to 
perform the necessary calculations to describe a blast wave of energy E q 
at one epoch, characterized by y $ ^. At that epoch, the program provides 
R s , t, y( r ) , x(r), T(r), T g (r), all ionic concentrations as functions of 
r, the column densities of all ions from near the explosion center to the 
edge, the total spectrum (distributed In 5 eV bins), and the total 
luminosity (emlssi vities use the calculated ion concentrations with 
the emission rate coefficients of Raymond and Smith, 1977). It 
furthermore calculates the spectral distribution of surface brightness of 
3 rays through the object, one passing through r * 0, and ones tangent to 
r * 0.5 and 0.9. Finally, for comparison, it recalculates the ionic 
concentrations, column densities, spectrum, luminosity, and surface 
brightness spectra which would result if the ions were in colllsional 
equilibrium with T e . The program can be switched to run either with 
coulomb heating of electrons, or with T e * T^ * T, but Is adiabatic in 
either case. As mentioned previously, the results should be appropriate 
for x $ > 2, y s > 3.5 so long as radiative cooling causes only negligible 
perturbations of the structure. For strong shock waves the latter 
condition requires t < t-j as given by equation (45). This condition Is 
modified by significant ambient pressure but can be roughly summarized by 
the requirement thatiSt must be small compared to E q . 

Eleven diagnostic runs were made with this program and some of the 
results are discussed in the remainder of this paper. 
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C. Runs Performed and Their Luminosity Evolution 

We are particularly Interested In applying this program to a study 
of the soft x-ray background, and this has biased our Initial set of 
calculations toward very low ambient densities and large scales. The 
actual runs which have been performed thus far are shown as 3 groups In 
Table 1. In group A, we followed the time evolution of an explosion of 5 
x 10®® ergs Into n Q ■ 5 x 10 cm , T Q ■ 5 x 10 K assuming that only 
coulomb collisions heat the electrons. This explosion was studied at 
four times bracketing t . In addition. It was restudied at two of the 
times with the alternative model, T ft » T^ « T. 

In group B, we studied three situations which differed initially 
only In their ambient temperatures. All had the same energy and ambient 
density and were examined at the same shock radius. All had t ~ t e ^ as 
defined for P 0 * 0 evolutions by equation (25). The differences in T 0 , 
however, resulted In differences in T $ , t, and all of the other 
parameters to some extent, which, as we shall see, caused variations of 
equilibration across this group. All three runs were made assuming that 
only coulomb collisions heated the electrons, but one was repeated with 
T ft ■ Tj a T. Within this group, run 31E was not entirely well-behaved, 
showi-sg some jumpiness In its properties. We incluJe the results here 
berause they contain information not present in other runs, but the 
detailed results should be considered as tentative. 

Group C consisted of only one run, with parameters similar to 
Idealized symmetric models of the Cygnus Loop (exclusive of the dense 
filaments). For both this example and run 4E of group A the equilibrium 
cooling rate was expected to be appropriate. As we shall see, this 
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prediction Is not supported In detail by the model results. 

Also shown in table 1, are the post shock temperature, the 
electron temperature for the Innermost zone, the calculated total 
luminosity, the total luminosity which would have been found If the Ion 
concentrations had been In colllslonal equilibrium with T e , and hand 
calculated estimates of the latter t*o quantities. 

The estimated luminosities were derived in the following way. The 
true luminosity is 


2 - 


/ 

o 


'4-rr it* 0 2 l] dR 


(56) 


whore L is the local cooling coefficient which depends on T e and all of 
the ionic concentrations. We assumed that this could be approximated by 

R. o o 

Z = L(T ' ) / s 4uR n dR 
e t 

(57) 

: 2.06 [ -j- nR s 3 n 0 2 ] L|T e ' ) 


where the compaction parameter of 2.06 was found from a numerical 
integration over Kahn's approximation to the Sedov solution. Here L ( T e * ) 
is the cooling coefficient at a representative electron . ?mperature. We 
chose T e ' as the post shock temperature or the innermost zone electron 
temperature, whichever was lower. 

To estimate the non-equilibrium 1 jmi nosit ies, we took L(T e ') from 
equation (50) with the result that for t < t^, 
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whereas the approximation to the RCS cooling function presented in 8IID 
was used for the collisional equilibrium luminosities to give 
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In comparing these estimates with the machine results, three points 
should be borne in mind. At early times (t < t e( ^), T e ' overestimates the 
representative electron temperature for the coulomb heating cases. In 
addition the assumption of a constant compaction parameter breaks down 
when T s /T q is not large, causing overestimated luminosities at late 
times. Thirdly, the very uncertain normalization of equation (50) was 
adjusted once the machine results were in hand to make equation (58) as 
accurate as possible. 

The estimates, by and large, are very good. No CE estimate errs by 
move than 40%; these estimates improve as t t e ^ as expected from 
above. For the T 0 * cases, the CE errors are less than 1 0%. for 
seven out of the 11 runs performed, the estimated nonequilibrium 
luminosities are also within 10% of the machine results, although it is 
clear from the run series 3E, 2E, IE, 4E, that where is a gradual 
Increase in the ratio of true to estimated luminosity. This drift could 
be removed if the typical value of E/I were taken as 0.37 rather than 0.5 
in deriving equation (50). This would then It^d to Z. a R^n^ '^(T^)^'^ 
t ^ When normalized to run 2E , this revision removes the drift in 
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group A, but alters the estimates for groups B and C by less than 5 and 
20 percent respectively, far less than would be required to overcome the 
significant discrepancies in 3 of the 5 runs. 

In addition, such an adjustment in the estimate leads to a false 
comfort that the luminosity evolution is fully understood in terms of the 
previous discussion. For run 4E, the compaction parameter is less than 
2.06 (as evivi- ,ced bv low T s /T 0 and the estimated exceeding the 
machine result}, and the luminosity estimate should be too high. For 
this run, however, the true luminosity is anomalously high just as it 
should be approaching the equilibrium rate. We believe that the reason 

for this depends on a detail: at the T s of run 4E, carbon would be fully 

ionized in equilibrium, but, in the flow, retains one or two K-shell 
electrons. The additional cooling due to these ions keeps the 
non-equilibrium rate higher than the equilibrium value until T s drops 
enough that carbon also has K-shell electrons in equilibrium. This same 
effect is seen exaggerated in run 21E where T s is high enough that oxygen 
performs this same function. Further discussion of this effect is found 
In Gronenschild and Mewe (1979). 

Finally the differences between the machine results and estimates 
of X for runs 3 1 E and 41E can be understood in terms of the effects of T 0 

on the ionization just inside the shock and is discussed briefly in the 

footnotes to the table. For cases 3 1 E and 21E, both with low T 0 , the 
magnitude of discrepancy may partly be due to poor resolution of the low 
stages of ionization near the edge, or to incomplete stripping of 
helium. 
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We have emphasized the discrepancies in this discussion, because 
they draw our attention to additional considerations In the physics. 

This would not be possible, however, if the overall quality of the 
estimates were not so high. The worst disagreement appears to be for run 
31E but, as we shall see, the agreement even for this run Is quite good 
if we limit our attention to x-ray energies. 

0. Some Details of Two Runs 

Both Itoh (1979) and Gronenschlld and Mewe (1979) have discussed the 
dynamic, thermal, ionization, and spectral structure of blast waves into 
negligible ambient pressure. We now present the details of two of our 
runs, reviewing tne principal conclusions of the earlier papers in the 
process. 

The runs chosen for this comparison are numbers HE and 11 from 
group B, differing only in their assumptions about electron heating. (As 
will be seen in § IV , these two runs have parameters compatable with the 
soft X-ray oackground model.) The dynamical structures of these two are 
identical and are shown in Figure 2a, b. The post shock compression 
factor is x $ » 3.05, the pressure enhancement y s ^ * 11.76. The effects 
of the ambient pressure on the structure are beginning to be felt at this 
epoch (R $ * 102 pc, R c • 106 pc). The mean temperature distribution is 
also the same for both runs, shown as a solid line in Figure 2c. In run 
11, this is also the electron temperature; but for run HE, the electron 
temperature is shown as a dashed line. Clearly this case is on the verge of 
equilibration of T , the interior electron temperature just slightly exceeding 
T s . From equation (40) this thermal structure is predicted to occur for 
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the p a 0 solution when t ^ 0.6 t whereas this run actually has t ^ 
o eq 

t as qiven by equation (25). Thus equilibration has been delayed by 
eq 

the lower compression and higher T $ brought about by the higher external 
pressure. Notice in Table 1 that run 31E, identical except for its low 
T 0 , has in fact equilibrated at this stage of development. 

The distributions of ionization stages for oxygen are shown in 
Figure 3a, b. The initial concentrations, equilibrium at 6 x lo\, were 

approximately 90% 0 , 7.5% 0 , and ?.5% 0 . The rapid disappearance 

+4 +5 

of 0 behind the shock is evident IN both cases; 0 follows a similar 

but slower decline. Had T 0 been lower; much of the oxygen In the dense 

+5 +6 

edge region would have been 0 , with 0 being the ionization stage in 

which the element then became stuck (in the sense of §110). For run HE, 
this sticking continues far into the interior, whereas for run 11, there 
is gradual conversion to the hydrogen-like ion, and even some complete 
stripping. This difference follows from the fact that the interior 
electron temperatures are much higher in the latter case, and were 
significantly higher yet In the past when this inner material was in the 
dense region near the edge. For run 11E, T g is nearly uniform In 
position and constant in time so that the ionization level depends 
primarily on /r^dt of each gas element. This integral is small near the 
edge, for newly shocked material, and small in the center for material 
which was shocked early but spent little time at high dens Hies. Both 
runs, in fact, show a peaking of the ionization level some distance 
Interior to the shock, but for run 11, this occurs farther inward than 
for 1 IE because of the higher electron temperature history of the former. 
The distribution of the equilibrium ionization structure of oxygen 
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is not shown, but it consists of a larger proportion of the hydrogen-like 
and bare ions. For fairly light elements at these temperatures, the 
equilibrium Is largely bare nuclei whereas the non-equilibrium case gets 
stuck trying to ionize the K-shell. For elements which are significantly 
heavier than oxygen, the K-shell is not reached in either case and the 
equilibrium and non-equilibrium distribution peaks can differ by several 
stages of ionization. For more complete discussions see Itoh (1979) and 
Gronenschild and Mewe (1979). 

Four of the surface brightness spect. jduced by these two runs 
are shown in Figure 4. These specific spectra are from a line integral 
through the center 


su,) 



. ( Ei) 

4tt 


dR 


(60) 


where c(E^) is the local emissivity of the energy bin characterized by 
E,j. These emissivities include line and continuum contributions summed 
in 500 bins of 5eV width. Thus a bin which stands a factor of 100 above 
the neighboring continuum has an equivalent width of 500eV. The four 
figures are the non-equilibrium (i.e. actual) and collisional equilibrium 
spectra of both runs. Note the differences in vertical scales. 

The major features of these diagrams can be summarized as follows: 

1) Lines dominate the luminosity for all four. 

2) The T 0 » T^ * T case (#11) has significantly more emission at 
high energies from the significantly higher electron temperature, and an 
upward curvature of the continuum from the hotter interior. 

3) The equilibrium continua show more pronounced recombination 
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edges from the high stages of ionization whereas the non-equilibrium 
continua show more pronounced two photon contributions from colli sional 
excitation of, for example, 0 + ® to a metastable level. 

4) The equilibrium spectra show significantly stronger lines from 
higher stages of ionization {appearing at high energies) than do the 
non-equilibrium spectra. 

As before, the reader is referred to Itoh (1979) and Gronenschild and 
Mewe (1979) for more complete discussions. 
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IV. Application to Model the Soft X-ray Background 

A. The Implied Parameters of the Local Emission Region 

The suggestion has been made, (e.g. McKee and Ostriker, 1977 ), 
that much of the soft X-ray background might arise from our being 
inside a very large supernova blast wave propagating in the hot, low 
density component of the ISM. This is a refinement of the notion of Cox 
and Smith (1974) that the background arises in typical material of the 
hot component, and this refinement lowers the inferred temperature and 
pressure requirements of typical locations at the expense of putting the 
solar system in an atypical region. Such atypical regions certainly 
exist, however, and the likelihood of being within one is not 
negligible. We have thus undertaken this project to explore that 
suggestion quantitatively. We stress, however, at the outset that the 
models we present do not correspond to those envisioned by McKee and 
Ostriker (1977) because the medium modeled does not have the large number 
of small clouds undergoing thermal evaporation as in their view. 

Although we dispute their point of view elsewhere for the ISM as a whole 
(Cox, 1979), for the present purposes we simply note that there is no 
evidence for there being such clouds, or even a typical amount of neutral 
hydrogen, within 100 pc from the sun. We thus model the explosion as an 
adiabatic blast wave. 

The gross features required for the disturbance come either directly 
from the observations, or from comparisons with the equilibrium 
emissivities of Raymond and Smith (19/7) for a hot plasma. Detailed 
discussions are found, for example, in Burstein, Borken, Kraushaar, and 
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Sanders (1976); Sanders, Kraushaar, Nousek, and Fried (1977); Nousek, Fried, 
Sanders, and Kraushaar (1981 ); Fried, Nousek, Sanders, and Kraushaar (1981); 
Cox (1977); Kraushaar (1979); and McCammon, Kraushaar, Sanders, and Burrows 
(1980). The basic results are that the temperature of the emitting gas (from 
the observed spectrum) must be ^10* K; the radial scale of the emitting 
region must extend to at least 100 pc (from the observed anticorrelation 
with moderate col umn densities of HI); and the likely emission measure 

2 -3 

/n^ dr is ^ 1 to 2 x 10 cm pc. Thus to first order, the density is 
-3 -3 -3 

about 4 x 10 cm , p/k ^ 8000 cm K. Finally, the thermal energy 

50 

contained in the nearest 100 pc is roughly 3 x 10 ergs, which is 
approximately the energy of one supernova explosion. It is this last 
result, of course, which makes this project tempting. 

A little manipulation leads to an estimate of the likelihood of 
being found within such a region. If the post shock temperature is 
roughly 10® K, the shock speed is about 300 km s V The age of an 

5 

adiabatic blastwave is t * (2/5)(R s /v $ ) or 1.4 x 10 years for R $ = 100 

pc, T « 10® K. If SN occur in the low density phase with a rate of 1 

per 100 years, there should be about 1400 in the Galaxy as young as the 

one envisioned here, most of them nearly as old and large as this one. 

65 3 

Thus such remnants occupy a total of about 10 cm of the galactic disk, 
or about 1 part in 40 of a disk with radius 15 kpc and thickness 200 pc. 
Their typical center to center separation, if positioned randomly in such 
a disk, would be about 500 pc. Thus the chances of a random point being 
located within such a disturbance are only of order 1 in 40 although the 
chance of being within 200 to 300 pc of the center of one or more is 
essentially unity. Finally, we note that the supernova reheating of a 
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preexisting cavity of this size bounded by a dense shell can resemble the 
model under discussion but have a significantly longer lifetime. Such 
models have been proposed for Loop I (Borken and Iwan, 1977) and for a 
region extending from Barnard's Loop to Eridanus (Reynolds and Ogden, 

1979) and may eventually be found to describe best our local region as 
well. 

B. Comparison of X-ray Observations witn Model Results 

The diffuse background data of the Wisconsion soft X-ray astronomy 
group is presented as count rates in a set of bands defined by a 
combination of window transparencies, counter response functions, and 
pulse height limits. Although detailed pulse height spectra are also 
obtained, the spectral resolution is such that most of the information Is 
carried In the band count rates. Three of these bands, arranged in order 
of increasing mean energy have been named B (for boron filter), C (for 
carbon filter), and M (for medium energy). For these bands, count rate 
maps of the entire sky are approaching completion. Tentative maps for 
80% of the sky can be found, for example, in McCammon et.al. (1980). More 
precise definitions of the bands are in Burstein, et. al. (1976). 

The B, C, and M bands are sensitive, respectively, to photons in 
approximately the ranges 100 to 188 eV, 150 to 284 eV, and 450 to 
lOOOeV. The count rates in these bands vary over the sky, but typically 
have values of order 50, 200, and 100 counts per second, respectively, in 
directions away from bright sources. The minimum values for these rates 
are about 30, 80, and 50 counts per second. 
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The soft X-ray astronomy group has programs which model the response 
of their instrument to input spectra, calculating the band count rates 
which would be seen for an Incident surface brightness spectrum. 

Bursteln, et. al . (1976), for example, have presented the results of such 
calculations using collisional equilibrium input spectra from Raymond and 
Smith (1977). The results are conveniently presented in band fraction 
diagrams (analogous to color-color plots or chromaticlty diagrams) and as 
the required emission measure to achieve a given count rate. 

The spectral results, summarized briefly, are that for temperatures 
below about 10^ K very little M band is produced but B and C bands are 
copious, the ratio of B to C increasing with decreasing T. For 
temperatures above about 10® K, M band and C band are bright, M/C 
increasing with increasing temperature, but little B band is produced. A 

temperature can be found at which B/M is as observed, but the predicted C 

band is then much too bright. Thus the equilibrium models are driven to 
the assumption that two temperatures are required, one higher than 10 6 K 
to produce most of the M but only part of the C, and one lower than 10® K 
to produce most of the B band and the remaining portion of the C rate. 

Nousek et al . (1981) have suggested that the hot component could be 
distantly located, such as in the galactic corona. From spectral studies of 
the lower energy bands In directions of varying hydrogen column density, how- 
ever, it appears that most of the emission sensed by the B and C bands arises 

20 -2 

closer than the nearest 1 x 10 atoms cm of hydrogen, generally within 100 

to 200 pc. The possibility that the M band emission alone arises from more 

distant regions Is consistent with the rather close correlation between 
variations in the B and C band sky maps, and the lack of correlation of 
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these two with the more uniform M band map. Thus, for the present 

program, our goal is to account for the B and C band rates, with the 

observed M band rate as an upper limit. 

The central chord surface brightness spectra calculated for the runs 

in Table 1 were divided by two and folded through the counter response of 

the soft X-ray counter to simulate the appearance of such blast waves 

viewed from within. The fact that we ar: unlikely to be found dead 

center is of little consequence since most of the emission arises near 

the edge, and surface brightness is distance independent. The resulting 

+5 

band count rates are shown in Table 2 along with the integrated 0 
column densities. Note that the count rates are shown for both the 
actual time dependent ionization structures and for the hypothetical 
equilibrium ionization distribution, for comparison. 

The total surface brightness of these models should be proportional 

2 5 

to«£/R . In fact, for all of the models with T Q = 5 x 10 K, the caroon 

band count •-ates of Table 2 are given within 6 percent by 

C = 15.6 [- g r — ] ( JOO P . c ) 2 counts s' 1 (61) 

10 35 erg s" 1 R s 

This accuracy must depend on the carbon-band-sensitive photons accounting 
for a large proportion of the total luminosity, or at least a nearly 
constant proportion. 

The exceptions to the above behavior are easily understood. For runs 
with low T Q (31E, 21E), the total luminosity is enhanced by the presence 
of low stages of ionization which do not radiate in the carbon band. For 
run 31E, however, equation (61 ) applies if the predicted luminosity from 
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equation (58) is used rather than the very high machine result. For run 
21E, the carbon band is intermediate in brightness between the values 
obtained from the predicted and machine calculated luminosities. 

Apparently the carbon band Is enhanced somewhat by the low T Q , but not by 
as large a factor as X as a whole. For the one run (41E) with a higher 
T , the converse behavior is found. The carbon band is less suppressed 
by the higher T Q than is£. Nevertheless, for runs 41E and 21E, 
equation (61) using predicted luminosities from equation (58) gives 
carbon band count rates within 14 and 37 percent respectively of the 
calculated values. Hence the carbon band rates in Table 2 are a reliable 
measure of the total surface brightness of the models with which to compare 
the B and M band rates. 

The degree of ionization equilibration has a strong effect on the 
count rates. This can be seen by comparing the CE rates with the actual 
count rates. In all cases, the non-equilibrium rates exceed the CE 
values, the difference being significantly greater for B band than for C 
or M. By following the run sequences (3E, 2E, IE, 4E) and (2,1), it Is 
clear that this non-equilibration enhancement decreases with time, as 
expected. Within the set (31E, HE, 41E). the enhancements are greater 
for larger T Q , but this can largely be accounted for just by the 
gradation in equilibration found earlier, with little additional effect 
depending on T Q alone. 

Since the C band rate evolves approximately as expected from the 
luminosity behavior, the ratios of band count rates M/C and B/C are 
sensitive to primarily two quantities: the temperature and the 

differential effect of equilibration. When the C and M band rates suffer 
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similar nonequilibration enhancements, as in the sequence (3E, 2E, IE, 

4E) the ratio of M/C should decrease in time due to the falling electron 
temperature. The B/C ratio, however, has competing tendencies. It 
increases as T fi declines, but decreases with Increasing equilibration. 

These effects are all seen in Table 3 which shows T e ', B/C, and M/C 
for all runs performed. For the run sequence (3E, 2E, IE, 4E), M/C 
declines monotonically by a factor of about 3 as y falls, while B/C 
shows a slight net increase, with fluctuations. For the runs (2,1) with 
T e * Ti» M/C drops by a factor of 3 while B/C increases by only a factor 
of 1.4. 

The B/C ratio thus shows little variation for the most part. For 8 
out of the 11 runs, it lies between 0.34 and 0.46. It is somewhat lower 
for the only 2 runs (2,21E) which have T ' > 2.4 x 10^ K and T g ~ T^. 

The only run with high B/C is 3 IE which has the lowest T e ', low T , and a 
somewhat questionable resolution of the edge In the model. 

The M/C ratio, on the other hand, shows considerable variation In 
the models; but this variation is of little interest simply because the 
ratio is always much lower than that observed in the actual data. For 
only 2 runs does M/C exceed even 0.10. Run 2 with M/C e 0.13 has both a 
high T e ' and T g ■ T^. Run 21E, although not designed to represent the 
background, never-the-less reemphasizes the required parameters for 
achieving high M/C. This run, like run 2, has high T ' and effectively 
T e =T.j . Furthermore, it is closer to ionization equilibrium than run 2 
which has a very high nonequilibrium enhancement of the C band. Thus 21E 
has a higher M/C than run 2 even though its T ' is somewhat lower. 

We are now prepared to compare the models with the data in a gross 
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sense. The first 6 runs of Table 2 all have C band count rates of about 
160 per second, typical of the sky data, Independent of their great 
differences. Similarly, the next 4 runs have C band rates of about 60 
per second, similar to the minimum sky rates, and independent of their 
differences. We have seen that these C band rates can be reasonably 
predicted by combining equations (58) and (61). For coulomb heating 
models which have not yet equilibrated their electron and ion 
temperatures, the resulting predicted carbon band rate can be written 


,7/5 


C ~ 116 ( t-nV" 

0.004 cm 0 10° K 


s ,1/10 




o 

30 " 


1/5 


(62) 


5 x 10 ergs 


where the Sedov solution (P 0 * 0) was assumed to remove the dependence on 
t and R , and T g ' was taken as the plateau value of equation (39). 
Similarly, for T e ; T^, the carbon band rate is approximately 


C 'v 86 ( 


"o ,4/3 , T s ,1/3 , E o 

it) ( — g — 1 ( — — eg — ‘ — 

0.004 cm J 10° K 5 x 10~ u ergs 


1/6 


(63 ) 


From these two equations, the Insensitivity of the carbon band rate to 
all parameters other than density is apparent, so long as T s ^ 10^ K is 
chosen as a prerequisite for fitting the spectrum, and p Q is not too 
large. Thus the models confirm the equi librium-calculated* estimate that 
n Q % 0.004 cm Is required of any model which attempts to reproduce the 
background surface brightness. We find here, however, that this is 
essentially independent of the radius of the blast wave disturbance, 
particularly for models with only coulomb heating of electrons. 
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In making a detailed comparison with the observational data. It Is 
necessary to consider the potential contribution from an extrapolation of 
the extragalactlc component which dominates at higher photon energies. A 

20 

power law extrapolation to low energies, assuming an Intervening N u * 4x10 

-2 n 
cm as typical, would contribute about 1 count per second in the B band, 

7 counts per second In the C band, and 21 counts per second to the M band. 

Only the M band Is altered appreciably and even there, most of the rate 

remains unexplained. The alterations which these contributions make In the 

band ratios are shown as B'/C' and M'/C' In Table 3. 

The observed B/C ratios range between 0.16 and 0.5, with 0.28 to 

0.42 being most common at low intensities (100 counts per second in C 

band) and 0.22 to 0.39 most common at higher intensities (200 counts per 

second In C band). Since the models have essentially this same range, 

the B/C ratio is only a weak discriminant for conditions in the blast 

wave. It is clear, however, from run 31E that unacceptably high B/C 

ratios result from low T e ' and low T Q . The ratio is expected to rise 

rapidly for T e ' < 10^ K, excluding these low temperatures from 

consideration. Conversely, a very low B/C ratio was achieved only by 

having T e ' > 2 x 10^ K, T e s T^ and progress toward ionization 

equilibrium. The ’atter condition, however, would require a density 

considerably higher than allowed by the overall surface brightness. 

In short, the observed B/C ratio is consistent with those models 

with the correct C band rate so long as T e ' lies in the range 1.3 to 2.5 
6 

x 10 K, whether or not the electrons have equilibrated with the 1o r 'S, 
and this consistancy extend to even higner values of T e ' if the 
electrons are heated only by coulomb collisions. 
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Finally, none of the models produce M band photons at a rate even 
close to what is observed on the sky. The simplest interpretation of 
this (other than that the models are wrong or inappropriate) is that the 
M band derives primarily from more distant regions, perhaps the galactic 
corona as mentioned previously. There is, however, an alternative 
explanation lurking in the results. We have seen that high production of 
M band requires both high temperature and approach to ionization 
equilibration, which in turn requires great age at the densities 
consistent with the C band rate. High temperatures and great age, 
however, are inconsistent with the models with ^ 100 pc since the age 
is roughly 0.4 R $ divided by the sound speed. 

On the other hand, both Loop I and the Barnard's Loop through 
Eridanus region mentioned previously appear significantly enhanced in M 
band, and both have been interpreted as the supernova reheating of a 
previously existing cavity. This interpretation depends in part on the 
apparent enclosure of these regions by HI shells which would not have 
been generated by the inferred supernova explosions at their observed 
states of evolution. Two points follow immediately: 1) although none of 

the models run here have produced a significant amount of M uand, there 
appear to be similar regions in the sky that do (so that we dare not 
assert that the background M band rate surely comes from more distant 
emission); and 2) the existence of a massive HI shell around a region 
allows confinement of the explosion to a finite region, possibly 
increasing the age enough that ionization equilibration is approached 
more fully. This possibility could be explored quantitatively by 
following the evolution of a blast wave into a region with a steep 
gradient in ambient density. We have applied the techniques of §1 and § 1 1 
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of this paper to this problem, but have not yet calculated the spectra of 
such regions. Since M band is preferentially produced in hotter regions, 
it may be necessary to include thermal conduction and evaporation into 
the interior in order to suitably model this component. 

C. The OVI Column Densities 

The OVI absorption measurements discussed by Jenkins ( 1978a, b) serve 
as Indicators of the general Interstellar environment as well as a 
constraint on the solar environment. With regard to the former, it is 

4-5 

possible that 0 generated by the blast waves in a very diffuse 

interstellar environment produces much of the observed absorption (e.g. 

McKee and Ostriker, 1977). We focus here however on the local 

constraint, noting that several stars with distances of order 100 to 500 

13 -2 

pc have N(OVI) £ 1.3 x 10 cm , implying that any blast wave with R s ^ 

100 pc in which the sun finds itself immersed can have N(0VI) no greater 

13 -2 

than about 10 cm . 

All but two of the runs performed have N(OVI) between 0.18 and 0.86 
13 -2 

x 10 cm in apparent consistency with this constraint (see Table 2). 

This agreement, however, was achieved by foresight. In all of these nine 

+6 

runs, the external ambient temperature was so high that 0 was the 

dominant Ionization stage before the shock even reached the material. 

1 3 -2 

Both runs with low values of T 0 had N(QVI) 'v 4 x 10 cm , distinctly at 
odds with the constraint. 

+5 

These machine results can be easily understood. The 0 ionization 
coefficient between 1 and 3 x 10^ K is roughly 10 ^ ( T/10 ^ K) cm^ s ^ 
(Summers, 1974). Thus an oxygen ion which suddenly finds itself in a 
high temperature region (temperature T $ for T * 


T^ cases) in which the 
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electron density is n^ % 4n Q will successively proceed through its 
various stages of ionization, spending a time 


t(0 +5 ) 


<o.v>n 
I e 


. 10 10 cm* 3 ,10 6 Kv 

^ {--y ) 


4rT 


(64 ) 
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as 0 . If this heating Is caused by a strong shockwave moving with 

speed v s , then the number of oxygen ions (with abundance A(0) relative to 
hydrogen) swept in per cm s is A (0) n v , a fraction f(< ♦ 5) of which 

U j 

are initially in the +5 or lower stages of ionization. The steady state 
column density of 0 ions behind the front is thus f {< + 5) A(0) n Q v s 
t(0 +5 ). For A(0) ^ 7 x 10 and v $ % 3 x 10^ (T s /10^ cm s \ the 

resulting column density should be 


N(OVI) 'v 5 x 10 13 ( 1 y-) 1/2 f(< +5) cm" 2 , 


(65) 


independent of n Q and the shock radius, and roughly in accord with the 
machine results, particularly for the T e = cases. The dependence on 
f(< + 5) is confirmed since for 21E and 31E, f - 1 whereas for all others 
except 41E, f * 0.1. For 41E, T q * 10^ K, f is even smaller and N(0VI) 
has its lowest value. In addition, the weak increase of column density 
with decreasing T s is apparent, as is about a factor of 2 higher N(0VI) 
for all T g < T^ cases. The latter is due to having a lower ionizaton 
coefficient in these runs than is appropriate for T , Finally we note 
that in all cases the column densities are significantly in excess of the 
equilibrium values, the excess decreasing with time. 

The consequences for the models are rather harsh. In order to hold 
the 0VI column density to an acceptable level, one of the following seems 



59 


to be required: 

a) The oxygen abundance Is about a factor of four lower than stated 

b) T q Is sufficiently high that most of the oxygen is in the +6 
stage; 

c) The ambient oxygen is found mostly in the +6 stage from past 
history even though T Q is now low; 

d) The post-shock temperature must be significantly higher to raise 

the 0 + '> ionization coefficient (which, however, has a maximum of only 8 x 

10~^cm2 s"^ at T % 10^ K); 

+5 

e) The 0 ions experience considerable photoi om zation; 

f) The rate of incorporation of new material proceeds at 
considerably less than the sound speed for the temperature to which they 
are heated; or 

g) The OVI constraint has been exaggerated. 

We have the following comments on these possibilities; 

a) The presence of oxygen-bearing dust could contribute to this 

effect, but the dust would have to survive in the post shock environment. 

5 

b) Any temperature greater than about 3 x 10 K is sufficient to 

13 -2 

make N(OVI) within the shocked region less than 10 cm . In this case, 

+5 

however the explosion serves to deplete a region of 0 , not to produce 

it. Thus there is the additional constraint that if the ambient medium 

-3 

(with density n Q * 0.004 cm ) occupies a significant portion of the 

+ 5 ^ Q _ n 

insterstel lar volume, it must have n(0 ) 5 2 x 10 cm. With A(0) % 7 

x 10"^, this implies that less than It of the oxygen is in the +5 stage, 

•> 

requiring T q : 8 x 10 K . At a just sliyhtly higher temperature, the 
ambient medium becomes a profuse radiator of soft X-rays and we revert to 
the model of Cox and Smith (1974) with no need for the local explosion. 
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c) This suggestion seems to be untenable. The radiative 

recombination timescale of 0 +6 with OgMx 10~ 3 cm" 3 and T e < 10 6 K Is 
7 

less than 1.3 x 10 years, but since only lit can have recombined without 

5 

violating the constraint on the ambient medium, less than 10 years can 
have elapsed since the temperature fell. This timescale Is much shorter 

than the approximately 10® years between strong blast waves in typical 
regions making recombination hysteresis an unacceptable hypothesis. This 
could be alleviated somewhat if the oxygen were left in the +7 or +8 
stage, but the implied high ionization stages of other elements would 
probably make production of B band photons difficult when the shock 
arrived. 

d) This may help somewhat in conjuction with either of 
possibilities f) or g). 

+5 

e) The photoionization timescale for 0 in the calculated 

15 +3 

radiation environments is of order 10 seconds, about a factor of 10 
too long to be of significance. 

f) This may be possible for models evolving in a pre-existing 
cavity with a strong density gradient and will be explored In future work. 

g) The constraint is only about a factor of four lower than the OVI 
colinn densities achieved for the low T 0 models. Since 0V1 is more 
difficult to measure for high temperature components which are broader 
and shallower than those actually seen, high temperature shocks could 
conceivably give profiles consistent with the measurements. 
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V. Summary and Conclusions 

In setting out to study the possibility that the soft X-ray 
background might derive from our being located inside a very large-scale 
blast wave, we soon realized that the external pressure was important to 
the details because the blast radii were in many cases comparable to the 
characteristic radius defined by equation (5). In §1, we managed to 
extend the approximation scheme of Kahn (1975) to find analytical 
expressions for the evolving dynamics. These descriptions are expected 
to be reasonably accurate throughout the structure for y > 6, R $ < 1 .27 
R ; and for the bulk of the material, which is located near the outer 

V 

edge, until at least y % 2, R % 2.5 R . 

3 ■'s * s c 

In §11, we discus' H the potential effects of electron heating and 
thermal conduction, settling on 2 adiabatic models which would reasonably 
describe most of the possible situations, particularly at late times. We 
also applied Itoh's (1978) formulation of coulomb heating to find 
analytic expressions for the electron temperature distribution of the 
evolving structure. These expressions are valid in detail only if this 
is the sole mechanism for electron heating and if thermal conduction Is 
negligible. It was shown, however, that with only coulomb collisions 
heating the electrons, thermal conduction would always be modest. The 
effects could, however, be just as important at t ^ t^ as are the non- 

negllgible external oressure effects at R s ^ R c * 

We were concerned that the blast waves we were studying might be 
beyond the point that radiative cooling upsets adi abat icity . Hence we 
restudied the question of the cooling timescale of a strong blast wave 
with particular attention to whether the equilibrium cooling coefficient 
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could reasonably be used In the evaluation. Our conclusion was that It could 
be used v*1th confidence even though colli si onal equilibrium of the Ions would 
not always be fully established. 

For strong blast waves, we were able to conclude that collision-related 
phenomena occurred In a fixed sequence: thermal conduction fades In Import- 

ance just as electrons are heated by coulomb collisions; later the colllslonal 
equilibrium cooling coefficient becomes a reasonable approximation, and 
finally significant cooling takes place. We expect this sequence to remain 
even when the ambient pressure is not negligible and conclude that we can 
tell whether radiative cooling is important in cases of Interest by observing 
the location of the results of detailed computer runs within this sequence. 
None of the cases In this study were found in this way to have left the 
adiabatic regime. 

In § 1 1 1 , we presented the method with which we were able to follow 
the history of an individual gas parcel through the structure so that the 
ionization equations could be integrated. We then outlined the computer 
program which used the methods of these 3 sections to model the blast 
waves in our study. For a set of diagnostic runs, the luminosity 
evolution was compared with simple models. It was found that the 
luminosity could be predicted to within a few percent in many cases; and 
when the estimates were found to be less accurate, the discrepancies 
could be rationalized. We presented the detailed structures and spectra 
of 2 runs to complete the picture before going on to our specific 
application to the soft X-ray background. 

In §IV, we presented the appearance of our diagnostic runs to the 
Wisconsin soft X-ray rocket payload, as count rates in the B, C, and M 
bands. For the C band rates, it was found that the rate for a particular 
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explosion varied very little with time and radius, and tha f this could be 
understood from the estimated luminosity evolution. The results for this 
band (with modest p Q ) are summarized in equations (62) for t < t^ and 
(63) for t > t or . We conclude immediately that the ambient 

_ 3 

density must be n Q * 0.004 cm to give the observed carbon band rate, 
although an implicit assumption, so that the dominant cooling falls in 
the C band, was that T s ^ 10® K. 

Once n Q has been determined by the C band surface brightness, the B 
band rate is still available to determine the evolutionary state (T s , R $ , 
t) for a given type of model. As it turns out, however, the ratio of B 
rate to C rate (B/C) is rather insensitive to the desired parameters, 
since it tends to increase as T e falls but decrease as the evolution 
approaches equilibrium. We are, however, encouraged by the fact that 
most calculated examples have B/C within the observed range. For the 
value of n Q determined by the C band rate, it appears to be possible to 
obtain too low a B band rate, only if electrons are equilibrated with the 
ion temperature by processes other than coulomb collisions and we look 
very early (when T e Is high). 

In contrast, the M band count rates produced by the models were 
never as large as the observed rates. We noted, however, that two large 
emission regions present in the soft X-ray maps, regions which could be 
more distant examples of such blast waves, are quite bright in M band. 
Thus, although it is possible simply to assume that for the general 
background the M band derives from some other source (such as the 
galactic corona) we would at the same time have to conclude that our 
bubble differs from the two others observed. In addition, the two 
others, at least, are not well described by the present models. 
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Neither of the other two bubbles, however, is expected to be 
described by the models. Both seem to have external hydrogen shells and 
have been Interpreted as blast waves evolving in pre-existing cavities, 
probably with strong density gradients. We believe that the presence of 
such a gradient acts in just such a way that M band enhancement will 
result and have begun to develop the models to test this hypothesis. 

The question still remains, even assuming the success of the next 
generation of models in producing M-band-sensiti ve photons, whether the 
general background derives its M band rate from our local bubble. This 
question probably assumes its greatest significance when we wonder what 
fraction of the count rate might be attributable to the galactic corona, 
and may translate to whether or not we are within a pre-existing cavity 

and whether the blast wave is interacting with a strong density gradient near 
the cavity boundary. We have only two comments: First, we may be located in 
a hole In the HI distribution that extends to about 100 pc or even further 
in some directions (e.g. Bohlin, Savage, and Drake, 1978). Second, we recall 
the initial reason for inferring that the local bubble had a radius as great 
as 100 pc: The observed anticorrelation between the count rates in both 

B and C bands and the HI column density along the line of sight was most easily 
interpreted as a displacement effect (Sanders, et al . , 1977; Cox, 1977). 

That Is, the hydrogen seemed to be intruding into the emitting volume, lowering 
the emission measure for soft X-rays. The clear implication, if this inter- 
pretation survives. Is that the blast wave is indeed interacting with strong 
density gradients In at least some directions. 

No strong conclusions, however, should be drawn from these comments. 
Although we seem to be in a cavity and our local bubble seems to be interacting 
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with the walls. It Is the B and C band rates which vary with this Interaction. 
The M band rate Is much more smoothly distributed. On the other hand, the 
M band rate Is also distributed more smoothly than the C band rate within 
Loop I (McCammon, et al. 1980). 

A relateJ problem not addressed by these models Is whether the presence 
of HI regions intruding into the boundary of the blast wave will in fact 
lower the B and C band rates, as desired, without seriously altering their 
ratio. 

Finally, we addressed a constraint placed on the blast wave ( and the 

ambient medium Into which it propagates) by OVI column density 

measurements. We found that for blast waves propagating into homogeneous 

cold material, the OVI co 1 umn densities exceeded the local constraint by 

a factor of about four. This could be remedied by having the ambient 

medium already highly ionized (mostly as 0 + ^), but then the ambient 

+5 

medium, alone, contained too much 0 unless the ambient temperature 
exceeded 8 x 10 6 K. With an ambient temperature just slightly higher 
than this lower limit, the blast wave is not needed to produce the X-rays 
and we revert to the initial model of Cox and Smith (1974). 

Other schemes for lowering the OVI column density were also 
discussed along with the possibility that the constraint was 
exaggerated. Two hopeful possibilities are related to the concept of 
there being a pre-existing cavity: for evolution in a strong density 

gradient, the OVI column density may be reduced to an acceptable level 
even for T 0 = 0; and, quite possibly, the "ambient medium" should not be 
thought of as extending throughout a large portion of the interstellar 
volume. It may, instead, be localized within large but discrete 
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cavities, most of which have been reheated sufficiently recently that 
+5 

little 0 is present except at the boundaries. 

Finally, we reiterate our successes; we have managed to model a 
blast wave in a homogeneous region of finite pressure, and to show that 
if the solar system were located within such a blast wave with n 0 % 0.004 
cm , E 0 ~ 5 x 10 50 ergs, and R s ^ 100 pc, the Wisconsin soft X-ray 
rocket payload would measure the B and C band count rates which it does. 

The explosion age would be just over 10 5 years. 
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TABLE 


LUMINOSITY EVOLUTION 


GROUP A: EVOLUTION IN TIME 

E g • 5 x IQ 50 ergs, n g « 5 x TO' 3 cm* 3 , T 0 • 5 x 10 5 K, R c « 86.0 pc, 
t^On 10 12 s, tj . 4 x 10 12 s, tj » 1.7 x 10 13 * 


Electrons heated only by coulomb collisions with Ions 


RUN 

♦ 

n 


HI 


T # (max) 

Do 6 *) 

if , (approx ) b 
(10 35 -ro S' 1 ) __ 

^ CF ,(*pprox) c 

3E 

55.7 

1.20 

4.82 


2.9* 

2.97 (3.4) 

1.19 

(.84) 

2E 

71.3 

2.20 

2.64* 


2.7 

5.27 (5.3) 

2.61 

0.8) 

IE 

83. 7 

3.18 

1.54 s 


2.4 

6.69 (6.2) 

4.40 

(3.9) 

4E 

102 

5.09 

1.34* 


2.6* 

10/ (8.6)® 

7.46 

(7.6)® 


Similar but with T e ■ T^ 

■ 

T 




2 

71.3 

2.20 

2.64* 



5.54 (5.3) 

1.67 

(1.8) 

1 

83.7 

3.18 

1.54* 


— 

6.74 (6.2) 

4.04 

(3.9) 




GROUP 

B: 

: DEPENDENCE ON T 0 





E 0 .5* 

10 50 ergs, 

% 

• 2.7 x 

10‘ 3 cm' 3 , R s - 102 pc, 





l eq 1 

4 x 10 12 s, 

*E 

: 6 x 10 

12 s, tj c 2.4 x 10 13 s 



RUN 

1 

m 

■ram 

m 

(io’ 2 >) 


T s 

(10 6 K) 

T e (max) JL(approx) b 
(10 6 K) (10 35 ero s' 1 ) 

<&{-.( approx ) c 
(IP 3 erg l_L 

31 E 

581 

3x1 0 3 

4.02 


1.27* 

2.2* 33.3 d (3.8) 

2.86 

(2.3) 

HE 

106 

5xl0 5 

3.85 


1 .93* 

2.3 3.97 (4.3) 

2.28 

(1 .8) 

41E 

83. 8 

lxlO 6 

3.72 


2.57 

2.4* 2.87 h (4.7) 

2.00 

(1.6) 



above : 

electrons 

heated only 

by coulomb collisions 





below: 

T . * T 1 * 

T 





11 

106 

5x1 0 5 

3.85 


1 .93* 

4.21 (4,3) 

Lii_ 

JUQ 


GROUP C: A DENSER EXAMPLE 

E q ■ 5 x 10^ ergi, n Q ■ 0.16 cm" 3 , T g ■ 1 x 10*K, R £ » 99.0 pc, 

t ^ 0.4 x 10 12 s, t r = 0.6 x 10 12 s, t, • 2.1 x 10 12 s 
e q t 1 

RUN R ? t T ? T^aiyx) JP, (approx) *CF ,(approx) c 

a (pc) (10 l2 sl (l6 6 K) (10 b K) (10 JS erg s' 1 ) (10 j:, erg s' 1 ) 

21E 21 0.60* 2.41* 6.5 e 179 1 (53 ) f 54.4 (51)' 



Footnotes to Table 1 


a) Electron temperature used In estimating luminosities 

b) From equation (58) 

c) From equation (59) 

35 -1 35 -1 

d) A substantial contribution, (29.0x10 erg s , versus 0.6x10 erg s 

for run HE, with hv<55eV) arises from low stages of Ionization which the 

program does not resolve well. The X-ray and EUV luminosity Is better 
represented by the approximate results. 

e) The plateau temperature Is not well defined for t > t a „. T a is still 

cl] C 

rising inward of our last zone. 

f) Example for which t « t£. 

g) Although nominally, t > t£, the external pressure has lowered the 
compaction parameter so that*t CE is over estimated. In addition. 
Incomplete ionization of the carbon K shell enhances the nonequi librium 
rate somewhat, delaying its approach to the equilibrium value. 

h) For large T„, two effects depress the cooling rate below the 
estimated result. One is the lower post shock compression factor; the 
other is the high preionization, resulting In higher (and less emissive) 
stages of Ionization in the cooling region. 

1) As In run 31 E, the total luminosity includes large, possibly erroneous 

35 1 

contributions at low photon energies (106x10 erg s’ for 5eV<hv<55eV and 
35 -1 

19x10 erg s for 55eV<hv<l 05eV) because the low T allows low transient 

o 

stages of ionization just Inside the shock. 



TABLE 2. BAND RATES AND N(OVI) 







































































































Table 3 


RATIOS OF BAND COUNT RATES 


RUN 

# 

(10 6 K) 

B/C 

M/C 

B'/C' 

M'/C 

3E 

2.9 

0.40 

0 06 

0.39 

0.19 

2E 

2.6 

0.39 

0.04 

0.38 

0.16 

IE 

1.5 

0.46 

0.02 

0.44 

0.16 

4E 

1.3 

0.45 

0.02 

0.44 

0.15 

2 

2.6 

0.28 

0.13 

0.27 

0.24 

| 1 

1.5 

0.38 

0.04 

0.37 

0.16 

31 E 

1.3 

0.54 

0.02 

0.48 

0.39 

HE 

1.9 

0.46 

0.03 

0.42 

0.35 

41 E 

in.. - — - — 

2.6 

— 

0.37 

0.03 

0.34 

0.33 

11 

1.9 

0.34 

0.07 

0.32 

0.33 

21 E 

D 

r - — ■ ■ 
0.21 

0.19 

0.21 

0.19 


Primed ratios Include 1 cps in B band, 7 cps In C band, and 
21 cps In M band as typical contribution of extrapolated and 
absorbed high energy extragalactic power law. 
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FIGURE CAPTIONS 


Figure 1. End result of adiabatic blast wave. Final density 

distribution, x^ * n fVn Q , and mass expelled from within R. 

Figure 2. Dynamical and thermal structures of runs 11 and HE. 

a) pressure; b) density; c) temperature, solid line Is 
mean temperature for both and electron temperature for 
run 11, whereas dashed line Is electron temperature for 
run HE. 

Figure 3. Distribution of Ionization stages of oxygen, a) run 11, 

b) run 11E. 

-23 

Figure 4. Surface brightness spectra in 5eV bins. Unit Is 10 erg 
-2 -1 -1 

cm s sr per bin. a) run 11, b) run 11 with collisional 
equilibrium Ion concentrations, c) run HE, d) run 11E with 
collisional equilibrium ion concentrations. 
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